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SYNOPSIS 


T. KUPPUSAMY, Ph.P., Indian Institute of Tedmology, 
Kanixir, September 1970, "Consolidation Characteristics of 
Non-homogene ous Clay layers under Time Pependent loading” . 

The present work involves the study of consolidation 
behaviour of soils subject to complex field conditions. Detailed 
analyses are presented for the theoretical consolidation 
characteristics of a thick layer and multi-layered systems 
subjected to constant as well as time-dependent loading conditions 
In this work a stratum with continuous variation of consolidation 
parameters with respect to depth is characterised as a thick 
layer while discontinuous variation of parameters with respect 
to depth is assumed in a layered system. The time- dependent 
construction type of loading is assumed as a linearly varying 
load during the construction period and as constant thereafter. 
Based on the theoretical investigations of the above field 
situations and loading environments, comprehensive conclusions 
on the consolidation behaviour of soils are presented. In some 
cases, simplified analyses using Terza^i’s original forinulation 
are carried out and compared with the more rigorous analyses 
presented here. 

t 

In a thick clay layer the effects "of the variations of 
the consolidation parameters are studied. The influence of the 
location of the drainage boundary on the corisolidation character- 
istics of a thick clay layer with single face drainage is 
investigated. An attempt is made to present a set of design 



charts for a thick clay layer for the assumed variation of 
parameters under time dependent construction type of loading. 

A new formulation is made for a non-homogene ous thick 
clay deposit consolidating with varying parameters during the 
process of consolidation. In this case the parameters are 
functions of depth and the pore pressure as consolidation 
progresses. The effects of the above general variation of 
parameters are analysed. This formulation leads to a non-linear 
equation, the solution of which is obtained by step by step 
numerical integration. 

In layered systems the effects of extreme ranges of 
parameters of individual layers are investigated. A set of 
design charts in a two layer system is attempted for a possible 
range of parameters. 

Rigorous mathematical analyses of layered systems under 
construction loading are presented. Closed form solutions for 
single and two layer system under construction loading are 
obtained. A three layer system is analysed by numerical method. 
The solutions are illustrated with examples. 

A simple experimental set-up is deviced for testing 
layered systems under construction type of loading. The experi- 
mental results of layered systems up to three layers under 
instantaneous as well as construction type of loading are 
compared with the results obtained by analytical procedures , 



CHAPIER I 


imoDucTioir to the erobeem 


Greneral 

Till today the essential basis of ssttlemenl; analysis 
is the one- dimensional consolidation theory as proposed by 
lerzaghi. (43) four decades ago. The general success of this 
theory in predicting the field settlements has been one of the 
strongest incentives in the earlier developments of the soil 
mechanics science. An adequate summary embodying all the further! 
advances in the consolidation theory is not possible in view 
of the extensive consolidation literature. However, appropriate ■ 

developments related to the present research work undertaken | 

! 

! 

are reviewed in subsequent chapters. Thou^ the state-of-art 

! 

of s.dttl^ment prediction has been streamlined by simplified 
theories, the heterogeneous nature of soil' deposits still offers i: 
considerable challenge in applying the idealised theoretical | 

j 

models to real engineering problems. At best, most of the | 

Ifebeories can serve as tools in arriving at a rational engineering 
judgement. With the developments in computation technique and | 
improved laboratory equipments, it has now become possible to 
develop soluticaas for complex cases. These solutions can offer 
better insight about the consolidation behaviour of soils under 
field conditions than those obtained from hi^ly simplified 
models. All the numerical computations in this study are done I 
with the help of the digital computer I.B.M, 7044. | 



CONSOIIDiffilON THEORIES 


The exfceneive research works carried out hy laany 
research y/orkers in the past following Terzaghi's theory can 
be categorised mainly as follows; 

a . I»!!at he mat i c al analy s is 

b. Testing procedure and material properties 

c. Application to field problems 

Ifeithematical analysis 

The overall process of consolidation consists of a 
fluid flow mechanism in a compressible material. From the 
continuity equation of compressible material, Terzaghi derived 
for one -dimensional case the basic differential equation for 
consolidation akin to the diffusion equation in heat-transfer 
literatirce as; 


G 




( 1 . 1 ) 


An extension of the above equation to three 
dimensional case leads to the equation, (27); 

■ Cy^u = ^ (1.2) 





where 
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Verigin (45), (46), solved three dimensional problems 
for different configurations of loading area using Eq. (i;2). 

The assumption of tJae constitutive Eq. ( 1 . 2 ) was shown to be 
not strictly valid by a more rational theory proposed by Biot 
(5). The right hand side of Eq. (1.2) is a volumetric strain 
rate rather than excess pore pressure rate. However, the above 
equation is strictly valid when reduced to one dimensional form. 
Biot ’s formulation is based on the assumption that a linearly 
elastic skeleton material is. filled up with pore fluid. Erom 
the basic equations of the theory of elasticity the following 
equations are fci“mqiated; 


+ 

i-2y^ 

aa 

„ a2i _ 

0 

(1.3) 

3x 

3X 

+ Q— 

1-2ft ■ 

33 

5y 

du 

0 

(1.4) 

+ 

1~2|a 

33 

dz 

^ _ 

” 5 z 

0 

(1.5) 


Erom the continuity equation of compressible material it is 
derived: 

9 0 

^ ( 1 . 6 ) 

where v^ , V 25 are displacements in x, y and z directions 
respectively, 

are Lame's elastic constants of soil skeleton, 
e is the volumetric strain, 
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C is the coefficient of consolidation, shd 

u is the pore pressure . 

The four unknowns , ^ 2 , and u are 
simultaneously solving the Eqs. 1.3 to 1.6 togeth®^ with the 
appropriate boundary conditions. The mathematical exercise 
is highly involved. Few plane strain problems by bhis method 
are solved, (6), (18), (19), (36). inalysis of situations 

by this method is hi^ly complex and almost imposs^^l®* 

In. the consolidation process the fluid flo'^ mechanism 
is often coupled with visco-elastic stress-strai^^"'*'^’^ relation- 
ship in order to explain the "Secondary consolidate®^' effect. 
The effective stress-strain-time relationship cah tie translated 
into rheological models composed of Maxwell and Kelvin units. 
Taylor (42) combined a linear spring and dashpot "t® explain 
the viscous effects in soils. Since then many investigators 
worked various linear and non-linear models to the 

time-def cr333ation behaviour of soils during the overall process 
of consolidation, (16 ), (20), (35), (47), (48). The rheological 
model analyst® can accommodate only the material properties. 

The lack of simulation of the geometric properti®® 
consolidation problems limits the scope, of such studies. 

Testing prqcQdure and material properties 

The e.ijjiple oedometer test plays an 
the evaluation of the consolidation parameters soils, 
hiversifiod results are obtained if the coh-vsntional 
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jesting procedure is not adhered to (11), (16), (40). The 
variations of the coefficients C, k and m during the consolidation 
process have been investigated by Anandakrishnan (2), Zeevaert 
'49). 7/hile recognising the limitations of the trends of 
variation of the parameters obtained in the laboratory it is 
relieved that the same trends will hold true quantitatively for 
rield situations. Therefore certain variations of these 
parameters are assumed in some of the problems analysed in this 
?esearch. However, further refined testing technique nay 
provide the true variations of the parameters. 

Lambe (15) and Davis (13) gave approximate methods of 
5omputing settlements to include three dimensional effect using 
iriaxial apparatus. By conducting drained triaxial compression 
iest for the appropriate stress path the volumetric and axial 
strains can be rioted. Using these values of strains the initial 
a.nd final settlements can be computed from the theory of elasti- 
city for axisymmetric as well as plane strain problems , Thou^ 
the above testing method is simple, it is extremely difficult 
bo successfully conduct the test without failing the sanples 
as they will be near the failure conditions, (21). 

Application to field problems 

In actual field situations, the consolidation parameters 
k, m and C vary spatially. Stratification and layered deposition 
are common features in the field. Even in homogeneous thick clay 
layers the variation of the parameters with depth can be 
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envisaged, (17), (34). However, in actual application to field 
problems the lack of appreciation of the manner of variation of 
parameters, the coefficients of permeability, c aapressibility 
and consolidation k, m and G respectively and the loading 
conditions resulted in resorting to oversimplified models. * 
Richart (28), Palmer and Brown (22), Sridharan and Hagaraj (39) 
have provided certain approximate method of analysing layered 
systems. The deviations that can result from the approximation 
are compared with the exact solutions in subsequent chapters. 

OBJECT I\m) SCOPS OE THE PRESENT IFVE'ST I CATION 

The scope of this research is to obtain solutions to 
several complex field conditions and loading’ environments so as to 
arrive at certain general conclusions on the consolidation behavi- 
our of soils subjected to these conditions and to compare with the 
results obtained by the simplified theory. From a review of 
several case studies involving layers with varying material 
properties it is seen that the rates of consolidation are evalu- 
ated with an almost total dependence on the coefficient of 
consolidation. One of the purposes of this study is to evaluate 
the effects of relative magnitudes of the variations in permeabi- 
lity, k and compressibility, m and the location of the drainage 
faces even though the coefficient of consolidation nBy still 
remain constant. A limited number of laboratory experiments 
simulating some of the above situations have been conducted and 
the results are compared with the analytical solutions. 
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Por th.e purpose of this study the following possible 
field situations are considered: 

G-eometrical configuration 

Most of the complex field situations can be idealised as 
sho\’m in Pigs. 1.1 to 1.3- The thickness ofthe compressible layers 
in these situations is small in comparison with the large loading, 
area at the surface. One dimensional analysis is valid for 
these cases. The discontinuous variation of” the soil properties 
with depth constitutes the layered system as shown in Pig. 1,2. 

A continuous variation of soil properties v^ith depth constitutes 
a thick layer as shown in Pig, 1»3. The variations are as shown 
in Pig. 1.4, 

lature of loading 

In the construction of buildings, embankments and other 
structures the loading always progresses gradually with respect 
to time. This construction type of loading is assumed as shown 
in Pig. 1.5. The total pressure, u^ is imposed in a time, t^ 
which is the construction period. The instantaneous loading is 
shown in Pig. 1.6. 

Variation of parameters during consolidation 

The variation of the parameters during consolidation is 

Ui 

assumed as function of (1 - — ) as shown in Pig. 1.7. 

o 
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Genera l variation of parameters 

In general, the parameters can vary with respect to 
depth and the pore pressure as the consolidation progresses. 

This is shown in Pig. 1.8. 

In Chapter II, an attempt is made to analyse the 
consolidation characteristics of a thick layer by numerical 
method. Different variations of soil parameters with respect 
to depth are investigated so as to arrive at comprehensive 
conclusions. The effects of relative situations of drainage 
boundaries and loading nature are also studied. The deviations 
in the results that may occur in such cases when simplified 
theory is used are clearly brou^t out. A new formulation of 
consolidation theory with general variation of parameters is 
presented. This forms a non-linear consolidation theory of a 
n on-homogeneous clay layer with variable coefficients. The 
solution of the non-linear equation is obtained by step by step 
numerical integration. Various regenerated cases can be derived 
from this general equation. The effects ofthe general variations 
of the parameters are analysed. If data regarding the variations 
of parameters in actual field conditions are available, the 
settlement analysis of such situations can be carried out by 
this present method only, to give reasonable estimates. 

In Chapter HI, layered systems under instantaneous 
loading are analysed. An attempt is made to study the effects 
of extreme ranges of parameters of individual layers. A set 
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of design charts in a two-layer system for a possible range of 
soil parameters is presented. The effects of the relative 
magnitudes of k and m- of the individual layers and the location 
of the drainage faces in layered systems are analysed. 

In Chapter IV, rigorous mathematical analyses of 
sin^e as well as two-layer systems are attempted for time 
dependent construction type of loading. The solution obtained 
for the single layer is identical with that of Schiff man’s ( 32 ), 
( 33 ), obtained in a different 7/ay. The aralytioal solutions for 
a t7/o~layer system under time dependent construction loading for 
single face as well as two face drainage are presented here for. 
the first time. The above solutions are checked with the numerical 
solutions of the above problems. Also the analysis of a three 
layer system is attempted by numerical method. The solutions" *- 

are illustrated with examples. 

In Chapter V,an investigation is taken for laboratory 
consolidation tests on layered systems for instantaneous as well 
as construction loading. A simple method of gradually applying 
the load to effect construction loading is device! . Tour types 
of soil are used to constitute the different layers in the 
tests. Samples with single, two and three layers are tested. 

The experimental results are compared with the theoretical 
solutions. 

Chapter VI summarises the whole investi^tion and the 
important conclusions are listed. 




FiG.l-? THICK LAYER FIG.i-4 VARIATION OF PARAMETERS 

WITH DEPTH 





CHAPTER II 


COISOLIDATION CE THICK CLAY LAYIE 
’■7ITH VARYING PAEAlETERS 

imODUCTION 

Demand for construction of structures on thick 
compressible soil is coumion in present-day engineering practice. 
The thick clay layer which is normally consolidating under 
its own overburden pressure ?/ill not have uniform soil properties 
throughout the depth.. The depth dependence of hydraulic 
properties as discussed by Dietrich (17) and Jackobson (4) shows 
that a linear ca? a non-linear variation of k with depth can 
exist in field. The variation of the consolidation parameters 
kj m and C in a thick clay layer can be due to the following 
factors : 

1. Decreasing void ratios at greater depths due to higher 
consolidation pressures. 

2. Diner materials which settle first at greater depths 
during sedimentation. 

3. Variation of soil itself with depth due to non- 
homogeneity. 

The variation of C with consolidation pressures are 
evident from laboratory test results, (41 )j (42). During the 
process of consolidation under an incremental load the void 
ratio changes. Therefore the parameters vary as functions 
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of pore pressure as consolidation progresses , (2), (31)? (34) > 
(49). It is also quite conceivable that even if the coefficient 
of consolidation, C is constant, different variations of k and 
m can exist. Hence an investigation of these factors will 
greatly enhance the understanding of consolidation behaviour 
of thick clay deposits. 

THICK CMY MYEE Y/ITH SPACE YARIATIOl OP k, m, AHD C 

Schiffman and Cibson (34) studied the consolidation 
of non-homogene ous thick layers for two way drainage only. In 
view of the variations of the parameters the effects due to 
one-way drainage cannot be extrapolated from the two-way 
drainage study. The location of the drainage boundaries will 
exert considerable influence on consolidation characteristics. 
The analysis conducted in this chapter for different variatioris 
of the parameters demonstrates the extent of approximation 
involved in using simplified approaches for such cases.. Also 
the extension of the analysis for construction loading is 
presented. 

Eormulation 

The general consolidation equation for a saturated 
clay mass in which there is a head generation is, (32); 

V K . V u + K V^u + Qm = m^ 


( 2,1 ) 



u 


where 


T7 - (— + — + — ) 
V - ^3x aY ^ dZ'*’ 


2 ^2 ^2 .2 

+ -L + -2^), 

dx dy az^ 


Q = the rate of head generation = , 


AP 


change in pressure, 


t = time, and 


u = excess pore pressure. 


For one dimensional case Eqn(2.l) reduces to 


.2. 


, 3"u dk da _ 

g^2 + az ■ 3z - ® ( at 


( ^ _ dC A P ) ^ 
dt 


( 2 . 2 ) 


Por constant k and instantaneous loading, Eq^. (2,2) simplifies 
to lerzaghi's original equation (which will be referred to as 
simplified equation hereafter): 


n . 2 _u _ 3QX 

t>„2 - at 

az 


(2.3) 


Pea: constant k and time dependent loading Eq. (2,2) takes the 
form; 


C ^ + R = M 

6z^ ^ 


(2.4) 
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in which R can be a function of space and time variables. 
Assuming the time dependent construction type of loading as 
shown in Pig. 1.6, the rate of imposition of excess pressure, 

R is equal to u^/t^ where u^ is the total imposed press-ure and 
t^ is the construction period. 

Por the case of thick layer with k and m varying 
with respect to depth alone the Eq. (2.2) is modified as: 

(2.5) 

Por the following variations of k and m, 

k = k^ (1 + rtx) and 
m = (“1 ± P^c) (2 .6 ) 

where 


X = the nbn-dimensionalised coordinate = 

H = the thickness of the layer and 
Z = the depth 


And using the following notations. 



C 


o 



T 



T 

o 


C t 
0 o 




and 


- 

c(z) - 


( 2 . 7 ) 
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from Eq. (2.5) it is obtained: 

(1 ±a + (±«) 9^" = 3P _ 1 

(1 +M) 9? (1+px) ix 


The boundary conditions are: 
1 . Drainage at both faces 


p(i,a?) 

P(0,T) 

P(x,0) 


0 

0 

0 


0 < T < T 

N O 


0 < X < 1 

N' V 


and 


2. Drainage at top face only 
P(0,T) =0 


9x 


0 < T < T 

V VO 


and 


(i,r) 


.= 0 


P(x,0) 


= 0 0 ^ X ^ 1 


After the construction is over the Eq. (2 
be modified as; 


(1 + ax) 

(1 + px) 


9^5* (+ ct ) 

“ = aF 


The boundary conditions are; 

1 . Drainage at both faces 

P^(x,T^) = P(x,T^), Solution of Eq. (2.8) 

P*(1,T) = 0 

^o v< ^ v< 


( 2 . 8 ) 


(2.9) 


( 2 . 10 ) 

.8) is to 

: T < CO 

( 2 . 1 , 1 ) 


P*(0,T) 


0 


( 2 . 12 ) 
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2. Drainage at top face only 

= P(x,T^)j Solution of Eq. (2.8) 

P^(0,'I') = 0 

;< T (2.13) 

3P^ I _ 0 

ax 1(1 ,T) - ° 


Por the case of instantaneous loading the following 
equation is to be used; 

(1+ax) a^p (+a) aP aP 

__Z + — = (2.14) 

("•iPs:) ('l+ps;) ax -aT 


The boundary conditions are; 

1 . Drainage at both faces 

P(0,T) = 0 

0 < T < o» and 

P(1,T) = 0 ^ ^ 

P(x,0) = 1.0 0 x^ 1 (2,15) 

2. Drainage at top face only 


P(0,T) 
ax , (1,T) 


P(x,0) 


0 

0 < T X and 

0 

1.0 0^ X ^ 1 


(2.16) 


The differences are evident by comparing the 
constitutive Eq. (2.14) and the simplified Eq. (2.3). In 
the Eq. (2.14), the coefficients of the differential terms 
are functions of the independent variable, x. The extra 
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second term in this equation reflects the Yariation of k. 

The Eq. (2.14) can be simplified for the following cases; 

Eor constant permeability, k (i.e., a = 0) 

-^ = (1 + px) II (2.17) 

ax 

Eor constant compressibility, m (i.e,,p= 0) 

(1 + ax) -2% +a f| = H (2.18) 

ax 

Einite difference equations 

Solutions of the Eq. (2.8) for various cases are 
complex by anal 3 rfcical approach. However, the equation is 
readily solvable with greater ease by using finite difference 
technique with the aid of a high speed digital computer. Ifeny 
numerical methods are available for solving the above type of 
equations, (29). The following finite difference method is 
adopted here. 

Using the finite difference mesh as shown in Eig. (2.1), 
Eq. (2.8) is written in the difference formas; 

(1 ±ax) (P 2 - 2 P^+P 4 ) (+a) ^^I'^o^ 1 

Ipx) ^^2 (1+px) AT T^ 


(2.19) 
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Rearranging the terms 




+ 1 1 - 


Ax^ 




AX 


T" 

o 


( 2 . 20 ) 


".'here 


A(x) 


I”* ± ax) 

TT+gx) 


B(x) = 


+ a 


(1 + Px) 


The boundary conditions are; 

1. At the drainage boundary, = 0 

2. At the impervious boundary, P 2 = P 


The stability and convergence of the above procedure 
depend upon the following conditions {34), 


AT 



< 

- dk 
dx 

< - ] 

- 2^t^ 


and 


• ( 2.21 ) 
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CASES STUDIED 

In order to study the effects of the variations of 
the parameters for different boundary as well as loading 
conditions, the following cases are analysed. 

Case 1 ; A tenfold increase in coefficient of permeability, 
k and coefficient of compressibility, m with constant coefficient 
of consolidation, C as shown in Eig. 2,2. 

Case 2; A hundredfold increase in coefficient of permeability, 
k and coefficient compressibility, m with constant coefficient 
of consolidation, C as shavn in Eig. 2.3. 

Case 32 A tenfold decrease in coefficient of permeability, k 
and coefficient of compressibility, m with constant coefficient 
of consolidation, C as shown in Eig. 2.4. 

Case 4! A. hundredfold decrease in coefficient of permeability, 
k and' coefficient of compressibility, m with constant coefficient 
of consolidation, C as shown in Eig. 2,5. 

In all the above four cases the coefficient of 
consolidation is constant. Analysis by simplified Eq. (2,3) 
will be a common practice for these situations. Therefore the 
study of the above eases will illustrate the differences that 
might exist between the results using the simplified approach 
and the present method. 
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Case 5* A tenfold increase in coefficient of permeability, k 
and a tenfold decrease in coefficient of compressibility, m 
effecting a hundredfold increase in the coefficient of consoli~ 
dation, C as shown in Fig. 2,6. 

Case 6; A tenfold decrease in coefficient of permeability k 
and a tenfold increase of coefficient of compressibility, m 
effecting a . hundredfold decrease in the coefficient of consoli- 
dation, C as shown in Fig. 2.7. 

The above two cases are analysed to compare with the 
analysis by the simplified Eq. (2,3) using average values of 
the coefficients. 

Some of the above cases are analysed, for the following 
conditions : 

1 . Drainage at both faces 

2 . Drainage at a single face 

3. Instantaneous loading 

4. Construction type of loading 

RESULTS 

Instantaneous loading 

The results of the computations for the cases studied 
are presented in the form of graphs in Figs. 2.9 to 2.16. The 
dotted curves shown in the figures are as per the simplified 
Eq. (2.3) whereas the other curves are as per the modified 
Eq. ( 2 . 14 ). The variations of the parameters are shown in 
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A clay layer draining on both faces and with constant 
coefficient of consolidation, C but with a ten or a hundredfold 
decrease of coefficients of permeability, k and compressibility, 
m, respectively show slower rate of consolidation by the analysis 
assuming space variation of the parameters k and m (Eqs2.14) as 
compared with the rate of consolidation by the simplified 
Eq, (2.3) as seen in Eig. 2.9. The tenfold decrease cons_olidates 
■faster than the hundredfold decrease. The pore pressure distri- 
bution along the depth at various times for the hundredfold 
decrease (Eig. 2.10) show a significant deviation from the 
simplified analysis. Therefore, even if the coefficient of 
consolidation is constant in a thick layer the estimation of 
the rate of settlement by the simplified equation will be in 
considerable error if the space variation of k and m is not 
considered. 

Eor a thick clay layer draining on both faces with a 
tenfold increase of k and a' tenfold decrease of m effecting 
a' hundredfold increase of C with depth, the consolidation - 
time curve is shown in Eig. 2,11. Also the case of tenfold 
decrease of k and a tenfold increase of m effecting a:, hundred- 
fold decrease of C with depth is shown in Eig. 2.12, Simplified 
analysis with average values of C at the mid depth and the 
^ithmetic average of the top and bottom values is done. The 
results of a similar analysis with a single face drainage 
located on the side with lower coefficient of consolidation 
and on the side of higher coefficient of consolidation are 
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shown in ligs. 2,13 and 2,14 respectively. The results by the 
simplified equation are considerably different from those 
obtained by considering the space variation of k and m. However, 
the simplified analysis with average C value at mid depth shows 
considerably less error than that by assuming arithmetic average. 

Another type of variation in k and m is visualised as 
a tenfold increase or a tenfold decrease of k and m but constant 
C with depth. A single drainage face is assumed to be located 
at the top only. In such a situation the drainage face is 
located on the side with lower values of k and m in the former 
case and on the side of higher values of k and m in the latter 
case. The results of the analysis of these two cases are shown 

in Hig, 2*15 along vvith the results by simplified equation, 

% 

These are designated as curves (l), (2) and (3) respectively in 
the Pig. 2.15. It is interesting to note that when the drainage 
face is located closer to the higher values of k and m the rate 
of consolidation is faster than if it is located on the side of 
lower values , Similarly a hundredfold increase of k and m 
(curve (4) in Pig, 2.14) sh 9 ws slower rate of consolidation than 
the above two cases. In these cases, it may be noted that the 
soil has higher compressibility farther away from the drainage 
face but lower permeability towards the drainage face. Thus the 
volume of water expelled due to larger volume change, farther 
from drainage face, meets with higher resistance to flow 
towards the drainage face due to the decreasing permeability. 
This causes the retardation in the rate of consolidation. 



Similarly, the accelerati-on in the rate of consolidation is 
effected in a reverse situation. Incidentally, the approximation 
involved hy simplified analysis is also brought out in this 
figure. Similarly for a thick clay layer with tenfold increase 
of k and m with depth and the single face drainage located on 
the higher values of k and m, the rate ofconsolidation . faster 
than if the drainage face is located on the side of .lower values 
of k and m (I’ig. 2,16). The simplified analysis for single face 
drainage does not bring out this difference. 

Construction loading 

The results of settlement analysis assuming construction 
type of loading for the various cases discussed under instantaneous 
loading are presented in T'igs . 2.17 to 2.29. As before the 
analysis is carried out for thick layer assuming the different 
variations of k and m with depth. For all these cases the 
consolidation-time curves for a set of T^ (=C^t^/fF/2) ) 
values corresponding to the construction periods, t^, are given 
in Figs. 2,17 to 2.25. These curves will serve as design charts 
for computation of settlements for the various situations shown. 
However, the range of T^ values covered is limited on account 
of enormous computation time required. For any intermediate 
construction periods within the chart, suitable interpolations 
can be made. Also in all these figures, the limiting conditions 
of instantaneous loading as well as the simplified analysis by 
Eq. 2.3 are shown in dotted lines for comparison. 
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In Pigs. 2,26 to 2.29 the curves are drawn relating 
the average pore pressure ratios, u/u^, where u is the integrated 
value of pore pressure throughout the depth, and the time factor 
T for various values of 1^. These curves show the pore pressiire 
build up and decay characteristics under construction loading. 
Here curve- (1 ) represents the simplified case wherein k, m and 
C are all constants throughout the depth. Curve (2) represents 
a tenfold increase while curve (3) represents a tenfold decrease 
of k and m with depth but having constant C. The drainage is 
allowed at the top only. The resistance . to flow causing 
retardation in the rate of pore pressure decay when the drainage 
face is located farther from higher compressibility but closer 
to lower permeability is observed here also as in the case 
of instantaneous loading. 

The simplified method shows different deviations depend- 
ing upon T^ values. The building up and decay characteristics 
of pore pressure depend upon the construction periods also. 
Realistic estimations of consolidation settlements for consti'uc- 
tion loading can be made only by the analysis with the modified 
equations presented here. 

THICK CLAY LAYER WITH GENERAL VARIATION CE k, m MU C 

Hon-linear consolidation probleuB have been attempted 
in the past with the assumption of linear e - log k and e~log P 
relations for homogeneous clay layer, (23), (24), (26). But 
the consolidation of non-homogene ous clay layer with variable 
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coefficients has not been attempted so far. In general, the 
parameters k, m and C can vary with depth and pore pressure 
as consolidation progresses. Some simple functions of k and 
m are to be assumed here in order to analyse the effects of 
the general variation of the parameters on the consolidation 
phenomenon. 

Pormulation 

The one dimensional consolidation equation for 
instantaneous loading for variable coefficients from Eq. (2.1) 
is ; 


a z 


( 2 . 22 ) 


when k = k (u,Z) and 
m = m (u,Z) 

Erom Eq. (2.22) it is obtained by chain rules 


K 


£ 

3Z 


» +1^ I 


5 + n 


8u 


au Bk ? 

az az 5 


Putting X = ^ I I* = ~ 

^o 


t* = 






we get 


(2.23) 


O ^ ”D ' ^ ^ Tr- *-1 •tT- _L 


as: 


2 ' k aP ''ax'' k ax ax k at' 


(2.24) 


Now the variation cf k and m is assumed as 


k = k^ (1 +ax) I 1-X(1-P)| and 


(2.25) 
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(2.26) 


(2.27) 





Substituting Eqs (2.27) and (2.28) into Sq. (2.24) and 
putting Oq = 5 ® 


We get ; 



lc^(1+ax)|1-^l1”!^) j dx 


_ (1+6x) I 1-H(1-P)) 
(l4ax) ] 1-5(l-P)? 
< •* 



(2.29) 


Simplifying we get? 


^ + A (»‘,?)('^)" + B 
3x‘“ ® 




(2.50) 


When e , 


A(x,P) = 'i_^(i_;p) 


= irfe ’ 
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C(a,p,if^ = 

The above Eq. (2.30) is non-linear. Closed form 
solution of the equation is not possible,. However, in the 
consolidation problems there always exists an open boundary. 
The values of P are specified only on the three boundaries as 
shown in Pig. 2.8. 

1 . At time t = 0 

2. At depth Z = 0 

3. At depth Z = H 

Taking advantage of the open boundary step-by-step 
numerical integration of non-linear consolidation equation is 
possible. However, very small intervals of integration should 
be chosen in order to get convergent results. Hence the 
computation time is large for solving such problems. 



Finite difference equations 

Using the finite difference mesh as shown in Fig. 2.1, 
Sq. (2.30) is written in difference form ast 



(2.31 ) 
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Therefore , 


P. = 


AT 


1 C(a,p,jr, lA ^ Pqj-) , „2 ■'■ 


P -2P +P P -P • 


iS. X 


P -P 

+ B(a,x) 4 ^ + P„ (2.32) 




The boundary conditions are: 

1 . At the drainage boundary = 0 

2. At the impervious bounday Pg = P^ 

The stability and convergence of the results depend 
upon the individual values of the parameters a , p , >!) |-i, Ax and 
iiT. However, the convergence of the particular cases studied 
is given in the Appendix. 

Regenerated cases 

The following regenerated cases can be derived from 
the general Eq. (2.29). 

1. Thick layer with space variation of k, m and C, which 
was discussed earlier . 

Here JT = |x = 0 

The Eq. (2.29) reduces to the form of Eq. (2.14). 

2 . A layer with varying k, m and C with respect to the 
pore pressure. 

Here a = p = 0 
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The Eq. (2,29) reduces to s 



+ 


'Jf 

i-y(i-P) 



qT 


(2.33) 


3. ^’Phen k alone varies as a function of pore pressure. 
Here jl- = 0 

The Eq. (2.33) reduces to; 


9^P 


8^ 


^ 1 _ 
2 ' i-)r(i-P) ~ i-xd-P) aT 


6P 


(2.34) 


4. When m alone varies as a function of pore pressure 
Here = 0 

The Eq. (2.33) reduces to; 


a^p 


ax' 


= 1 -^(1-P) n 


(2.35) 


Schiffman (32 ) gave an approximate niethod of solving an equation 
similar to Eq. (2.34). Scott (31 )» studied similar variation 
of Eq. (2.35). 


GASES STUDIED 

In order to understand the effect of the general 
variations of the parameters, the following cases are analysed. 

Case 1; A layer with variation of k alone during the process 
of eonsolidation. So here « = P = j-t = 0 
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Case 2t A layer with variation of m alone during the process 
of consolidation. So here a = p = =0 

The above two cases are taken for analysing the 
relative effects of the individual variation of k and m. 

Case 3’ A layer with k and m varying as the same functions of u. 
S o here a = p = 0 and ?" = ju. . 

k = k^ J 

Therefore , 



Case 4? A thick layer with k and m varying as the same functions 
of the depth and pore pressure. So here a = p and Y = p-. 

k = k^ (Hux) 

m i= (l -^x) I 1-p(l"P) ^ = m^(l-Kxx) I 1-'* (1-P) j 
Therefore, 

0 = t = ^2 = 0 

m 0 

In the above cases 3 and 4 the value of C is constants 
throughout the depth and dirring consolidation. Hence, the 
usage of simplified approach can be studied and the results 
can be compared. 
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Case A thick layer with general variation of k and m. 

So here a p V K 

A comparison of the results with the simplified 
analysis using average coefficient values will demonstrate the 
approximation involved. 

All the above cases are analysed for instantaneous 
loading and for two way drainage . 

RSSUITS 

The results of all the cases studied are presented in 
Rigs. 2.30 to 2.33. 

It is seen that the decrease in m during the process of 

consolidation assuming constant k, shows faster rate of pore 

pressiire dissipation than the decrease in k during the process 

of consolidation assuming constant m (Rig. 2.30). These effects 

K 

are obvious because the coefficient of consolidation, C = — , 
increases when m decreases and thus faster consolidation takes 
place. Similarly, C decreases v\fhen k decreases and hence 
slower rate of consolidation is seen. But the retardation of 
consolidation by a decrease in k is more than the acceleration 
of consolidation by the same decrease in m. This shows that 
the effect of the variation of k is more than that of m. 

Ror a. clay layer having the same coefficient of 
consolidation throughout the depth and during the consolidation 



process, the decrease of k and m result in slower rates of 
consolidation compared to the simplified analysis (Pigs. 2,31 
and 2.32). 

A general variation of k and m is effected hy taking 
the following values : 

a = - .80 ; p = 4.0 5 -5 = .25 and ^ = .50 

This shows that the permeability is decreased to one 

fifth and the compressibility is increased to five times the 

values of k , m with respect to depth. At the end of consoli- 
0 0 

dation process the permeability is reduced to three fourths and 
the compressibility is reduced to half the initial values, 
k(Z) and m(Z), The average value of C at the mid depth and 
at 50 percent consolidation is .233 0 ^. Pig. 2,33 shows that 
the use of simplified Eq. ( 2 . 3 ) with average value of coefficients 
results in considerable variation of rate of consolidation as 
compared to the general variable coefficient Eq. ( 2 . 29 ). 

CONCLUSIONS 

The following conclusions summarise the findings of 
this investigation: 

Por a thick clay layer even thou^ the coefficient of 
consolidation C is constant throu^out the depth, the rate of 
consolidation is very much affected by the individual variation 
of k and m with respect to depth. A tenfold and hundredfold 
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decrease in k and m result in slower rates of consolidation as 
compared to the simplified analysis. 

The location of the drainage boundary has significant 
influence on the consolidation process of thick clay layers 
with single face drainage. The consolidation is faster when 
the drainage boundary is located on the side of hi^er k and 
m values . 


In a thick layer under construction loading the pcxre 
pressure build up during the construction. and decay thereafter 
depend upon the construction periods also. Analysis of such 
cases with the modified Sq, (2.14) presented here will give 
better estimates of rates of consolidation. 

A set of design charts for a thick layer under 
construction loading are presented. If appropriate field 
situations could be idealised under the various cases discussed 
here the design charts will be useful for predicting the field 
settlements.. 

A new formulation taking into account the general 
variation of the consolidation parameters is presented. This 
forms a non-linear consolidation theory of non-homo geneous 
thick clay deposit with v^iation of parameters during the 
process of consolidation. Various regenerated cases can be 
derived from this general formulation. 
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The effect of individual decrease of k during consoli- 
dation is found to "be more than that of m. 

Even though the coefficient of consolidation is constant 
through-out the depth and during consolidation, the rate of settle- 
ment is influenced by the individual variations of k and m. 

In general, using average coefficients with the 
simplified Eq. (2,3) yields very approximate estimates as. 
compared with more refined analysis . 
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FIG-2-33 effect OF GENERAL VAPiA^lON k AND rr, 



CHAITEE III 


CONSOIIDiTION OP LjffERED SYSTEM 
UNDER INSTANTMEOUS lOADING 

INTRODUCTION 

The effect of several variables on the prediction of 
time rate of settlement using Terzaghi's original formulation 
has engaged the attention of many research workers in the past. 
However, the layered systems wh:|.ch are common occurrences in 
field have received only limitel attention. Several empirical 
approaches aimed at the use of transformed sections or averaging 
the soil properties yield approximate estimates, (22), (28), (39 )» 
The only closed form solution for a two-layer system was presented 
by Gferay (14). -Analytical approach for more than two layers is 
extremely difficult. However, the use of numerical technique 
with the help of high speed digital computers offers greater 
scope and facility for obtaining more dependable prediction of 
time rate of settlement for layered systems. Selected cases of 
layered systems. are analysed here for possible ranges of soil 
properties so that their effects on the consolidation phenomenon 
can be studied. 


FINITE DIFFERENCE METHOD 

Abbot (1 ) solved a three-layer system by finite 
difference method for a particular configuration of the layers. 
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Raymond (25) used finite difference solutions to coiapare with 
his experimental results of two- layer system, Barden (3) 
employed the numerical method to compare the test results for 
radial flow in layered systems. Davis (12) gave time factor 
contours at fifty percent consolidation in a two-layer system. 

An attempt is made here to present design charts for a two-layer 
system. Also, typical cases of three-layer system are analysed 
to demonstrate the errors involved in estimation by simplified 
approach. 


The simplified consolidation equations for layered ‘ 
system are: 


0 . 


a u. 

1 

a t 


(3«1) 


where, the subscript j denotes the layer number. 

The above equation is written in the finite difference 
form using the mesh as shown in Pig, 3*1 as follows: 


Therefore 


where 


u.-u C 

cl 


“1 = 


C. 6t 


(6Z.)< 


(3.2) 

(5.3) 

(3.4) 
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The boundary conditions are: 

(1 ) At the drainage boundary = 0 

( 2 ) At the impervious boundary U 2 = u^ 

( 3 ) At the interface of the two layers 

Consider the finite difference mesh on the interface 
boundary such that the points 0 , 3 and 1 fall on the boundary 
and points 2 and 4 on layers A and B respectively. 

In layer A: Prom Taylor's expansion we have: 






where OA denotes the derivatives at the point 0 towards layer A. 


Therefore 


(■H) = 

0Z^ OA 


(szr 


I U 2 -u^+ffl(||)^J (3.6) 


■The time derivative is approximated as : 


u.-u 
1 o 


'^>04 - ^ 


(3.7) 


Substituting the Eqs . ( 3 . 6 ) and (3.7) in the consoli- 
dation equation: 


^u 9 u , 




(3.8) 
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( 

We get — ~ < 
hir I 


Ug - Uq + f; 




(3.9) 


c; r7 / u H '' 


- "^0 . 

C^ 6t ^0 


(3-10) 


In layer Bs Similarly it can be derived in layer B. 


- 6Z(^) 


OB 


Ui-Uq 


(6Z)‘ 


+ ^0 " ^4 


(3.11) 


Now substituting the boundary conditions! 






(3.12) 


and ttLus equating the Eqs. (3.10) and (5.11) and simplifying, 


" 1=^0 "■ 


2>^ 

0 A 


/i '^A,. , ^A 

"o +^"2 


(3.13) 


where A = 6 t/(6 Z )' 


It can be readily checked that when k^=kg and 0^=0^, 
the Eq. (3.13) reduces to the original form of Ecj. (3.2). 


Alternatively this boundary condition can be derived 
by equating the rate of change of volume ahd flow in two half 



67 


sublayers, (25). Ihe expression by this procedure is obtained 
as % 

~ ^0 ^ 5 . 14 ) 

where 

" 3 , 3*1 = 2 = 3*1 “ 3*1 “" 3*/'“3 ^^ 3 ™ 3*1 «" 3 * 1 > 

and the subscripts denote the layers . 

The condition of Eq. (3.16) reduces to that of 
Eq. ( 3 . 13 ) for the particular case of 6 Zj^ = 6 Z^, 

The above finite difference method is employed in the 
explicit method to solve the simaltaneous difference equatiors 
obtained, noting the stability of this procedure depends upon 

p 

the criteria (C.6t/( 5 Z.) ^ 1 / 2 ). Here the computation is 
carried on in a partitioned manner. The time interval is 
divided into six parts and the space interval into ten parts. 
The values computed at end points of one panel is fed into 
the starting points of the next panel. As corsolidation 
progresses the mesh sizes are varied to save computation time. 

A general programme is written to vary all the parameters. 

Using the finite difference method, a fev/ selected 
cases of the two and three- layer systems as shown in Eig. 3.2 
to 3.5 are studied. The results of the computations of various 
cases are presented in the form of graphs in Figs. 3.6 to 3.35. 
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TWO- LAYER SYSTEMS 

Effect of the soil parameters k, m and 0 

In order to study tlie effects of the soil parameters, 
a two- layer system with layers of equal thickness is taken. 

The ratios k 2 /k^ , and varied. The pore pressure 

distribution along the depth at various time factors is shown 
in the figures . 

Layers with same compressibility 

The results show that in a two- layer system with same 
compressibility but with unequal permeability drained on both 
faces, the layer with higher permeability has faster dissipation 
of pore pressure (Eigs. 3.6 to 3.10). This behaviour can be 
visualised intutively. The volume of water expelled in both 
the layers is same because of the equal compressibility. In 
the layer with higher permeability, the water has to flow 
faster resulting in a higher rate of consolidation. However, 
when one-layer has relatively larger permeability, say, of the 
order of hundred times or more as that of the other layer, 
most of the consolidation in the earlier stages is essentially 
within the layer with hi^ier permeability. The rate of 
consolidation in this layer is very nearly the same as if it 
is a single layer with one way drainage. Therefore, if the 
rate of consolidation of a two-layer system with constant 
compressibility but with higher permeability ratios ( >100) 
is to be estimated in approximate steps it is suggested that 
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the initial rates be obtained simply for the layer with higher 
permeability for one-way drainage. The rates of settlement 
at longer periods can then be obta.ined by computing the rates 
for the layer with lower permeability treating it with two- 
way drainage, with its initial time taken at the end of the 
consolidation period for the layer with higher permeability. 

The above approximate method is clearly demonstrated 
by the example shown in T’ig, The total settlement of 

the two layers is assumed as 10 inches so that the settlement 
in each of the layers is 5 inches. The graph is drawn between 
the settlement and the time factor. 

Layers with same coefficient of consolidation 

Por a system with hvo layers of equal thickness and 
same coefficient of consolidation, with two face drainage 
there is no effect due to variation in k and m (Figs. 3.12 and 
3 . 13 ). This is because that the interface boundary condition 
is unaffected as 8u/ 3Z at the interface is zero and the flow 
is taking place from, the centre towards the top and bottom 
■ surfaces symmetrically. 

In as much as the coefficient of consolidation is the 
same in all the layers the normal tendency will be to treat the 
’ layer system as one single layer and use the simplified equation 
to obtain the rates of settlement. The fallacy in such an 
approach is demonstrated herein. Take for instance a two-layer 
system drained at the top face only with layers of equal thickness 



70 


and same coefficient of consolidation but the layer adjacent 
to the drainage boundary is having relatively lesser values 
of k and m. In such a situation the rate of consolidation in 
both the layers is found to be much slower than what is obtained 
by treating the system as a single layer (Pig. 3*14). In the 
figure the dotted c\irves are drawn as per simplified analysis 
assuming the layer system as a single layer with thickness 
equal to the total thickness of the layers. Ihe consolidation 
ratios at various time factors of the individual layers are 
obtained from the area of the pore pressure isochrones of 
the time factors in the corresponding thickness. It is noted 
here that the volume of water expelled due to larger volume 
change farther from the drainage face meets with higher 
resistance to flow tov/ards the drainage face through the layer 
with lesser permeability. Ihis causes the retardation in the 
rate of consolidation. 

The reverse situation occurs when the layer adjacent 
to the drainage boundary has relatively higher permeability 
and compressibility. The consolidation is found to be faster 
than that obtained by treating the system as single layer 
(Pig 3.15). 

Layers with same permeability 

In a two-layer system with a two-face drainage and 
layers of equal thickness and same permeability but with 
different compressibility, the layer with less compressibility 
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drains faster than the layer with higher compressihility 
(Mgs. 3.16 and 3.17). 

layers with different permeability and compressibility 

The combined layer action is more pronounced in a two- 
layer system with layers having different permeability and 
compressibility (Pigs. 3.18 and 3.19). Ihe pore pressure 
distribution shown in the figures demonstrates the combined 
influence of the permeability and compressibility ratios on the 
consolidation rates. 


Design charts 

Por a two-layer system draining on both faces, an 
attempt is made to obtain design charts for a possible range 
of soil parameters. The following non-dimensional parameters 
are chosen for the computations : 







The parameters cfand^are varied from 1 to 100. The 
consolidation ratios u^ and of the two layers are presented 
in Pigs. 3.20 to 3.31 at various time factors. Por values of 
"a and ^ from 1 to .01 , the position of the layers can be 
reversed and the same design charts can be used. Suitable 
interpolations can be made for any intermediate values of the 
parameters. 
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Illustration 

Por a case of a = 40 and |3 = 5 the consolidation ratios 
obtained from the charts are given in Table 3.1 


Table 3.1 

Consolidation of the two layers 



“■i 


T 

^2 

11.5 

.001 

30.0 

.08 

15.0 

.016 

43.0 

.16 

19.0 

.024 

47.0 

.24 

21 .5 

.032 

63.0 

.48 

81 .5 

.240 

84.0 

o 

00 

• 

97.0 

.480 

99.0 

1 .60 


The design charts for a limited range of the parameters 
are presented here. Since the computation time required is very 
large for extreme ratios of parameters, all the ranges are not 
attempted. 

THREE-IAYER SYSTEM 

layers with same coefficient of consolidation 

A three-layer system with two face drainage and layers 
of equal thickness and same coefficient of consolidation but 
with middle layer having relatively larger values of permeability 
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and compressibility is considered. In such a situation the 
rate of consolidation in al3. the three layers is found to be 
much slower than what is obtained by treating the system as 
single layer (Fig, 3.32). The drainage faces are located 
adjacent to the two outer layers having lower values of k and 
m than the middle layer. Therefore, the volume of water 
expected due to larger volume change farther from the drainage 
faces meets with hi^er resistance to flow towards the drainage 
faces throu^ the layers wdth lower permeability. This causes 
the retardation of flow. 

In a reverse situation where the top and bottom layers 
have values of k and m relatively larger than the middle layer 
the consolidation is faster (Fig, 3.33). 

The above cases are also studied for a single face 
drainage at the top. When the middle layer has relatively hi^er 
permeability and compressibility the rates -of consolidation are 
slower than the rates computed by simplified analysis assuming 
the system as a single layer (Fig. 3.34). In the reverse 
situation when the top bottom layers have relatively higher 
values of k and m than the middle layer, the rate of consolidation 
in the earlier stages is faster and in the latter stages slower 
than the fate predicted by simplified analysis assuming the 
system as a single layer (Fig, 3.35). The combined influence 
of the location of the drainage boundary and the relative values 
of k and m of the layers is observed here. 
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COUCLUSIOIS 

In a two-layer system with two face drainage and layers 
of same thickness and compressibility, if one of the layers has 
relatively very large permeability, most of the consolidation 
in the earlier stages is essentially in the layer with hi^er 
permeability. This phenomenon is utilized in suggesting an 
approximate method for computing the rates of settlement in such 
cases . 

A tv/o-layer system with two face drainage having layers 
of same thickness and coefficient of consolidation can, be treated 
as a single layer. However, if the drainage is on only one side 
it will be erroneouF. .to treat this system as single layer despite 
the same coefficient of consolidation. The rate of consolidation 
in such a case will be governed by the relative values of k and 
m of the two layers, and the location of the drainage boundary. 

Design charts for a limited range of soil parameters 
for a two-layer system with two way draixxage are presented. 

Dor a three-layer system with layers of same thickness 
and coefficient of consolidation whether draining on single or 
both faces it will be erroneous to resort to simplified analysis 
assuming the system as a single layer. In this case the rates 
of consolidation are governed by the relative values of permeabi- 
lity and compressibility of the individual layers and the location 
of the drainage boundary. 
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Using Terzaghi’s one -dimensional consolidation theciry 
various boundary value problems were successfully solved by 
many research workers in Soil Mechanics, The effect of some 
of the parameters ' involved in the consolidation phenomenon 
has been studied in the past . These researches were confined 
to instantaneous loading only. However, in the construction 
of buildings and embankments application of load requires 
considerable time and the loading progresses gradually. An 
approximate method of predicting the settlement of a single 
layer under construction loading was given by Terzaghi (41 ) 
based on intuitive arguments. Schiffman (33) obtained analytical 
solution for a single layer under construction loading. Here 
the mathematical susceptibility of single as well as two-layer 
problems is to be taken advantage of, for obtaining closed form 
solutions. However, for more than two layers the analytical 
approach is highly complex. Therefore, numerical technique 
involving either direct integration or application of convolution 
integral has to be employed for solving systems containing more 
than two layers. An attempt is made here to solve a three - 
layer system under construction loading by finite difference 


method . 
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SIlGIiE LAYER 


Analytical solution 

A single layer with two way drainage is taken as shown 
in Rig. 4.1 . The one-dimensional consolidation equation under 
time dependent loading iss 


C 



+ R = 


6 u 
5 t 


where Z is the space variable. 


The boundary conditions are: 


u(0, t ) = 0 

0 

< 

*s 

t 


u(2H,t) = 0 

0 

x< 

t 


o 

ti 

o 

0 

v< 

z 

^ 2H 


(4.1 ) 


(4.2) 


where 2H is the thickness of the layer and t^ is the construction 
period. Using the following transformations: 


P 




and 


d = 


EH^ 


n 


Eq. (4.1) can be transformed as | 


8~u , j ^ ('A 

0 p 


To solve the Eq. (4.3) with the boundary conditions corresponding 
to Eq. (4.2), the following method is used (38): 



t02 


Taking the equation — ^ ~Q^ (4.4) 

S p 

and using the teclinique of separation of variables, the solution 
to Eq. ( 4 . 4 ) which vanishes at p = 0 and p = 2 is; 

2 2 

00 ^ ni 

n It 4 

u = D^Sin(^^p)e (4.5) 

n=1 

where are constants and n are integer values. 

Now the solution to Eq. (4.3) is assumed to be of the form; 

00 

u = 5 : iZf(T) Sin ^ p • . (4.6) 

n=1 

where ^(T) is a function of T. 


TC 

Expressing d = I' E^Sin — ^ p 

n=1 


(4.7) 


where E^ are the Eourier Sine coefficients evaluated as 
(n = 1 , 3 ...). Substituting Eqs. (4.6) ajad (4..7) in 
Eq. ( 4 . 3 )» the following ordinary differential equation is 
obtained; 


djZ?(T) 

dT 


2 ^2 

£_>2L_ 0(T) 


4d 

n-jt 


(4.8) 
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Therefore : 


0{T) 


16d 




(1-e 


2 2 
n 7t 


4 


T 

), n = 


3*. (4.9) 


And the solution to Eq. (4.3) is 


u(p,T) 


16 u^ 
0 

tt^T 


- I 

o n=1 ,Sj . . 


•^(1-e 

XT 


2 2 
n % 


4 


T 

) Sin p 

(4.10) 


v/here 



Eor the post construction period t^ <: t , the consolidation 
equation is: 


^u^ _ r) u* 


(4.11) 


in which u* is the point pore pressure during the post- 
construction period and the boundary conditions are; 


u*(p,T ) = u(p 5 T ), Solution obtained in the previous 

case (as per Eq. 4.10) 


u*(0,T) = 0 

u*(2,T) = 0 


( 4 . 12 ) 
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The solution of Eq. (4.11) satisfying the boundary conditions 
of Eq. (4.12) is; 


16u 

u*(p,T) ^ ^ 


2 ^2 
n It 


T 


•d-e 


7t>l^ 1,3.. 


n' 


) Sin ^ p 


16u 


0 


T 


1 


Tt^T . , n' 

o 1,3.. 


(1-e 


2 2 

- ■7 ^- (T-T ) 
4 o . 


nit 


) Sin p 


(4.13) 


Defining the average pore pressure values u as the 
integrated value of the point pore pressure u(p,T) over the 
depth, the average pore pressure ratio ~ is given by; 


2 2 
n It 


^ K Uci.e 

u„ it^ . , „4 

0 o 1 ,3 . . 31 


), 0 <T <T^ (4.14) 


11 


32u 


0 o i ^3 • 


^(1-e 

n ■ 


2 2 
n^ Tt^^ 


T 


32u 


T It 
o 


i: 


•d-e 


2^2 

n_^(T..To) 


), T <T< 


1 ,3.. 


n 


(4.15) 
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'The solutions given by Eqs. (4.10) and (4.13) are 
identical with the solutions obtained by Schiffman, (33) 
using a different procedure. Also solutions obtained using 
finite difference approach yield identical results (Eigs. 4.4 
and 4.5). 

TWO LAYER SYSTEM WITH TWO-WAY DRAINAGE 
Analytical solution 

A tv/o~ layered system with two-way drainage as shown 
in Eig. 4.2 is taken. The one-dimensional consolidation 
equations for construction type of loading are; 


d^u. 

C. j + 

' az^ 

3^1 

(4.16) 


a^u^ 

cu^ 


Op — 1 + 


(4.17) 


in which the subscripts 1 and 2 denote the corresponding Soil 
layers 1 and 2. 


Using the following transformations; 
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the Eqs. (4.16) and (4.17) are transformed as; 




SP 


2 + di 


0^1 

TTF" 


(4.18) 


aS ^ . „2 

+ '*2 - I * sT 


(4.19) 


the boundary conditions are; 

u^(0,T) = 0 
u-| (p,0) = 0 
U2(p,0) = 0 

U2(1 + V ,T) = 0 


0 ^ T ^ Tq 

0 p .< 1 

1 ^ P 0 + 

0 ^ T ^ T 


(4.20) 


At the interface p = 1 


u 


^ (p,T) = U2(P5T) 


9u.j 

"lT~p 


= k, 


P=1'^ 


9u 2 
2 s p 


P=1 i 


The method of .solution is given in the Appendix. The solutions 
of the Eqs'. (4.18) and (4.19) readily satisfying the boundary 
conditions of Eq, (4.20) are; 

A^m 

” 2d., 1 , . 

u.^(p,T) = y Ej^(1-e ) Sin A^p Sin A^pV (4-2l) 

i=i -s 



107 


i=1 


-A?T 

1 


) Sin Sin (l + V-p) 


( 4 . 22 ) 


where A^s are the roots of the equAciont 


(cTcos A Sin Ap.V + Sin A cos A^iV ) = 0 (4.23) 


and = 


O' Sin A^uv + Sin A^ 
^AVA-Sin^ A^ + (TASin^A^ 


(4.24) 


Defining the average pore pressinre and U 2 as the 
integrated value of u^ and U 2 over the respective depths, the 
average pore pressure ratios are given by; 


u °° “4^ 

10 ^ 

= ^ ^D.(l-e 

% °i=i 


) Sin (1-^os A^) (4.25) 


-2 _ i. 

T. 


0 . . 

1=1 




'F- (1-e ) Sin A^ (1-cos 


(4.26) 


Dor the post -eonstruct ion period < T < » the consolidation 

equations are; 


(4.27) 



1-oe 


B^u* 


2 "^^2 


_ u. ^ 

~ r dT 


( 4 . 28 ) 


The b oxiridaxy conditions are 


u* (p, T^) = (p, Tq) 


Solution obtained in the 
previous case 


4(0, T) 


u*(l+V, T) 


U-o (Pj ) 


= U2(p, Tq) 


Solution obtained in the 
previous case 


At the interface p = 1 


u*(p, T) 


u*(p, T) 




p=lf 


-2 Tp ' 


(4.29) 


p=l4' 


The solutions of the equations (4.27) and (4.28) 
readily satisfying the boundary conditions given by Eq« (4.29) 


are : 


u*(p, T) = 


«> -A^T 

2d. i 

' 75^ ^ 

i=1 ^i 


2d- 

^ ~r 

i=t ^ 


-A7(T-T ) 

X ^ o 


) Sin A^p Sin (4.30) 
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2 dp 

u^Cp, T) = I: i’j_(1-e ) Sin Si]; A^(l+V~p.) 


- aJt 


.=1 


2 d . 


.=1 ^iK 


5 % ^1 ( 1 -« 


- AJ ( T - To ) 


) Sin A^ Sin A^^(1+ V -p) 

(4.31) 


The average pore pressure ratios are given by: 


lil _ 1 _ 

- m 

Uq 0 


j 3 5'i(1~e 
i=1 - 1 


- A^T 


) Sin A^p. (1-cos A^) 


^ y . 2 _ i .( i_e ) 


1 


0 . . A: 
1=1 1 


Sin Aj^ p (1-cos A^) (4.32) 


“ -A?T 

^ 91,9 ^ s 

= rfr* — F’j ('i-e ) Sin A. (l--cos A.U’v*) 

0 ° 1=1 1 


U 


0 


r. 3^ p. (n-e 

1=1 4^'> 


-A‘(T-Tq) 


) Sin A. (1-cos 

•it X* 


(4.35) 
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TWO-LAYER SYSTEM WITH ONE-WAY LRAT N AffE 
Analytical solution 

A two-layered system with single face drainage as 
shown in Eig. 4.3 is taken. Correspondingly the boundary 
conditions are; 

(p,0) = 0 


5u^ 


0 


!p=o 

U2(1+v,T) = 0 

U2(p,0) = 0 


(4.34) 


At the interface p = 1 

U.J (p,T) = U2(p,I) 




^1 ap 


= k, 


P=1 f 


2 ~^ 


p=1 ^ 


The solutions of the Egs . (4.18) and (4.19) readily 
satisfying the boundary conditions given by Eq. (4.34) are: 


2d. "'4'^ 

u^(p,T) = 2 E^(l-e ) cos A^p Sin 4^|lV (4*35) 

i=1 h 
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2d, 


2 


U 2 (p> 5 )) = V, ~2~2 ) *^08 Sin A^|Ji(l+V -p) 

i=:1 


(4.36) 


where A^s are the roots of the equation: 

(ffSin A Sin AjAV - Cos A Cos AjAV) = 0 
Cos A; 


and 'F. = 


(o-A. Sin^A^V+AuuCos^Aj) 


(4.37) 

(4.38) 


The average pore pressuxe ratios are given ty: 



1 

oo 

V 

__ 

^0 

II 

O 1 

/ 

i=1 



oo 


1 


Uq 

^ ^0 

z 

i=1 


-A^T 

P^(l-e ) S 


-A?T 

1 


in A^fJivSin A^ (4.39) 


^ ^ Z P. (1-e ) Cos A. (1-Cos A,l^v) (4.40) 




Por post-construction period Tq ^T , the solutions of 
Eqs, (4.27) and (4.28) readily satisfying the hotindary 
conditions similar to Eq. (4.29) are: 


2d 


A^T 

1 


u*(p,T) = >: E^d-e ) Cos Aj_p Sin 

i=1 \ 


°° 2d. ~A^(^-’^0^ 

y ) Cos A^p Sin A.^V 

i=1 h 


(4.41) 



1t2 


u, 


2d. 


-A^T) 


*(p,T) = 7 - 2 ' " 2 ) Cos Sin A^(l+v-“P)p. 

i=1 -^iF 


2d 


^ S', (1-e 


-Aj(T-To) 


^ 2 2 i 

1=1 


) Cos Sin A^^{1+ V ~p) 

(4.42) 


The average pore pressure ratios are given by; 

oo -A?T 

^110 ^ 

— ^ E ^ F. (1-e ) Sin A, Sin A.LIV 

-Lq P 1 1 11 

^•0 ^ i=1 1 


° 1=1 ^ 


Pei-io) 

e ) Sin A.^ Sin Aj_|XV (4.43) 


_* 


u 


2 1 2 


-A^T 


— ^ ^ F. (1-e ) Cos A. (1-cos A lAv) 

Un 0 . . Al’p ^ 1 11 

0 1=1 i« 


^ S Pi(l-6 

° i=1 




) Cos A^(l-cos A^aV) (4.44) 


Illustration 

The solutions for a tw-o-layer system with single face 
drainage are illustrated by the following problem. 

The soil parameters used for the computations are: 



16 


h^ 


= 5 aiKi 


k. 


= 8s 
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For these values Eqs . (4.39)j (4.40), (4.43) and (4.44) 
are evaluated using IBM 7044 digital computer and the results 
are presented in Figs. 4.6 to 4.9. 

Solutions obtained by analytical as well as finite 
difference methods for a two- layer system with two face drainage 
as given in Figs. 4.10 and 4.11, show very good agreement between 
the two methods . 


COMPARISON OF TERZAGHI'S I5ETH0D WITH RI&OROUS 
ANALYTICAB SOLUTION FOR SING-LE LAYER 

Terzaghi’s method is based on the intuitive argument 
that the settlement at the end of construction period is equal 
to the settlement which would have occurred at half of the 
construction period if the total load had been applied instant- 
aneously at the starting. Similar arguments can be extended to 
any time t and the settlement S^ is found as follows; 


C* _ Q 

"t - t^ ^t/2 


0 < t < t 


0 


(4.45) 


^t ’^(■t-'^0/2^ 


tQ < t < CO 


(4.46) 


The average pore pressure ratios of the analytical 
solution are converted to settlement ratios as follows; 


"^0 u 


0 f "fco 


(4.47) 


0 


= (■»-—) ^ ^ 


0 


(4.46) 


'0 
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The results of TerzagJii's approximate method and the 
rigorous analytical solution are compared in I'ig. 4.12. Por 
practical usage Terzaghi’s method proves to be a good approxi- 
mation of the analytical solutions. However, deviations occur 
for higher construction periods (Tq > .4). 

THREE-LAYER SYSTEM 
finite difference method 

The one-dimensional consolidation equation for layered 
system under construction loading is: 

a^u. au. 

C. — ^ + R = ^ (4.49) 

Where the subscript j denotes the layer number. There 
will be as many number of the above equations as there are layer 
The above equations can be transformed into the difference 
equations using the finite difference mesh as shown in Pig. 3.1 
as follows : 

ui 

where B. 

J 

The boundary conditions are: 

1 , At the drainage boundary = 0 

2. At the impervious bounday U 2 = 


u 

= (1-2B. ) Uq+B. (u2+u^) + ^t:^ 


(4.50) 


C. 3t 
J 

(SZ)^ 
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3 . At the interface of the layers • 
As given in Chapter III, 


u. 


^1 = ^^t^ (4.51) 


where = 2 (B .m.f.Z . )/(m. .^-Z f,Z ) 

^d,0+1 " ^^®D+l“o+l’'^D+1 “i+1 ^^d + 1^ 

Using the above finite difference method a three-layer 
system is solved for the following case (Eig. 3.4) s 


^2_1 

~ ’ C^ " 10 ^ IT - 3 


k, 

k 


h. 


- - - 2 f H - 5 


H = h^+h2+h3 ; T = Gt/(H/2)^ 


The average pore pressure ratios vs time factor is 
shown in Eig. 4.13 for a particular Tq = .50. Similar curves 
can be obtained for different time factors Tq corresponding 
to the construction periods tQ. However, it is to be noted 
here that computation time will be enormous if extreme ratios 
of the parameters of the layers are encountered. Because very 
small time intervals are to be taken to satisfy the stability 
criteria in order to get converging results . 
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CONCLUSIONS 

In evaluating the consolidation characteristics of 
soils, it is more realistic to assume time dependent cons- 
truction type of loading. The closed form solutions presented 
here for a single as well as tv/o- layer system will be helpful 
in predicting the time rate of settlement for field conditions. 
The computer involvement in evolving the final series of the 
analytical solutions is negligible compared to the time required 
for the finite difference method. When extreme ratios of soil 
parameters of the layers are encountered the computation time 
by the finite difference method is enormous . 

For systems of more than two layers the finite 
difference method is the only method v/hich can be easily 
employed for evaluating the consolidation characteristics, 

Terzaghi’s approximate method of estimating the time 
rate of settlement for a single layer under construction 
loading proves to be a good approximation of the rigox'ous 
solutions. 



IL 


■ z 

FIG.4.1 SINGLE LAYER WITH TWO WAY DRAINAGE 


C 2 k2 h2 

Cl ki hi z 


e 


FiG.4.2 TWO LAYERED SYSTEM WITH TWO WAY 
DRAINAGE 


C2 k2 h2 

C! k] hj 

y - y 

IG. 4-3 TWO LAYERED SYSTEM WITH ONE WAY 
DRAINAGE 






: \ 


. 4.4 


' n;. 


r- 










o 



I- 

CC 

g 

U- 

UJ 


I— 


O 


O 


O 


O 
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EXPEEII.!EmL SET-UP 

A consolidation apparatus where the loading is applied 
by lever system is used. Special 3 inches dia. and f inch high 
brass fixed type eoedometer rings, with grooves -cut in the 
rings, so that one can seat over the other exactly, as shown 
in Pig. 5.1, are used. The instantaneous loading is applied 
by placing weights at one end of the 3. ever arm. A gradual 
increase in loading at the end of the lever arm is accomplished 
by a continuous increase of vs'ater level in a plastic container 
at the end of the lever arm as shown in Pig. 5.2. A supply 
tank maintains constant head for uniform flow. The rate of 
flow is regulated by means of a control valve fixed in the 
supply tank. This arrangement typically represents a gradually 
applied construction type of loading in one-dimensional 
consolidation tests. 


soils TESTED 

Pour different soils are used as layers in the 
experiments. They are localy available Kanpur clay, commercial 
Pentonite, Kaolinite and a mixture of Kanpur and Bentonite clay 
called Tabe clay. The properties of these soils are given in 
Table 5.1. The grain size distribution curves are shown in 
Pig. 5.3. 
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Table 3»1 
Soil properties 


■Soil 

L L, 

P.I. 

S p . Gr . 

Gasagfando 

classification 

Kanpur clay 

44 

24 

2.60 

CL 

Bent onit e 

680 

632 

2.85 

GH 

Kaolinite 

53 

23 

2.62 

CL 

Table clay 

92 

64 

2,68 

CH 


TESTING mCCEDTIRE 

Tbe soil samples are prepared in the form of thick 
slurry. Thin filter papers are placed at the interface of the 
layers. To ensure complete saturation the soil samples are 
left with small back pressures for one day. To minimise the 
effects of side friction high graded silicon grease is applied 
on the sides of the oedometer rings. Two of the important 
limitations of the testing procedure are; 

1. the temperature effect, and 

2. the secondary effect 

The effects of temperature variation is minimised by 
conducting all tests in an air-conditioned room where the 
variation of temperature is between 18° and 22 °C. The secondary 
effect especially in tlie case of Bentonite and Tabcc clay are 
bound to be predominant. However, for the purpose of this study 
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the consolidation parameters obtained from the test results 
of a single layer is assumed to be applicable to the same soil 
when it is part of a malti-layer system. The interpretation 
of test results are based on this approximation 

EESUITS 

The consolidation characteristics of the soils used 
are found out by routine one -dimensional consolidation tests. 
The variation of the coci l'icient of consolidation y/ith stress 
levels are shown in Fig. 5.4. The void ratio- pres sure curves 
of the soils are given in Fig. 5.5. Using the coefficient of 
compressibility, a^, obtained from the slope of the void ratio 
pressure curve, at various stress levels, the coefficients of 
permeability are calculated and are presented in Table 5.2. 
These va^lues are subsequently used to obtain the analytical 
and numerical solutions of the various cases tested for the 
appropriate stress levels . 

Table 5.2 

Permeability values k cin/mln. 

level Table clay Kaolinite Kanpur clay Bentonite 
T/Sq . 


0.25 

7.50 

X 

10"^ 

6.10 X 

10“^ 

2.10 

X 

10“'^ 

2.00 

X 

10“^ 

0.50 

2.30 

X 

10-^ 

2.04 X 

10“^ 

0.72 

X 

10”"^ 

1 .00 

X 

10"^ 

1 .00 

0.90 

X 

10“^ 

0.44 X 

10”^ 

0.23 

X 

10”^ 

0.41 

X 

10“^ 

2,00 

0.36 

X 

10-^ 

0.17 X 

10"^ 

0.09 

X 

10“^ 

0.20 

X 

10“^ 

4.00 

0.25 

X 

10-^ 

0.07 X 

10"^ 

0.04 

X 

10""^ 

0.15 

X 

10~^ 
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In the case of layered systems under InstantLarLeous 
loading the experimental time-consolidation cuirves are compared 
with those obtained by approximate method as well as by more 
rigorous analysis outlined in Chapter III. The approximate 
method is based on transforming the thiclmess of individual 
layers (39) as follows; 





w^here H! 

D 

c. 


is the transformed thickness of jth layer, 

is the original thickness of 3 th layer, 

and C. are the coefficients of consolidation of 
J 

the first and jth layers respectively. 


n 


The total thickness H = 2 H. is assumed to consolidate as 

3 = 1 - 

single layer wath coefficient of consolidation C-j , 


CASES STUDIED 

Eor the various cases discussed in the previous 
chapters, certain idealised assumptions regarding the relative 
values of the parameters k, m and C were made with a view to 
study the effects of the variables either individually or jointly 
on the consolidation characteristics of layered systems.- While 
the analytical study is valuable to provide a basic understanding 
of the effects of the variation of the parameters and the use 
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of different interpretation technique, it is not possible to 
strictly duplicate the idealised extreme situations in laboratory. 
However, certain general trends and conclusions obtained from 
the analytical study can be meaningfully verified from the 
laboratory tests using the available soils. The following 
cases are taken up for the study: 

1 . Single layer under construction loading 

2 . Two-layer system under instantaneous . loading 

3. Two- layer system under construction loading 

4. Three-layer system under instantaneous loading 

5. Three- layer system under construction loading 

DISCUS SIOHS 

Construction loading 

Single layer: The test results of single layer with t\^/o face 
drainage under construction loading are given in Digs. 5.6 to 
5.14 hy dotted lines. The tests are conducted for various rate 
of loading as denoted by Tq, the time factor corresponding to 
the loading time t^. Theoretical solutions of these cases by 
the methods outlined in Chapter IV for the corresponding Tq 
values are shown by solid lines in these figures. 

Reasonable agreement of the test results with the 
theoretical prediction is seen (Digs. 5.6 to 5.10). However, 
noticeable deviation occur in Bentonite (Dig. 5.10). In this 
test the total imposed pressiire is 2.9 T/sq.ft. It is found 
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that the value of C under instantaneous loading itself varies 
fron 9.9 to 5.7 x 10"" ‘ cm /min. Therefore, taking these values 
of C, theoretical curves arc drawn as shovm in this figure. 

The test results are seen to fall within the limiting theoretical 
curves . 


In iPigs, 5.11 to 5.14 curves are drawn between the 
average pore pressiire ratios — (obtained from measured settle- 
ment ratios) and the time factor T at various Tq. These figures 
demonstrate the building up and decay characteristics of the 
pore pressure under construction loading. 

Tv/o and three-layer systems ? The results of the two and three 
layer systems under construction loading are shown in Figs. 

5.15 to 5.23. The analytical solutions of the various cases 
as given in Chapter IV are shown in all the figures for 
comparison. The laboratory experiments are conducted for 
various rates of loading. The results obtained are shown for ■ 

I 

the Tq values corresponding to the loading period tQ. 

The comparison presented here in all the figures shows i 
the amount of error involved in predicting the rate of settle- 
ment of layered systems under construction loading by theore- ! 
tical procedure. Hov/ever, the general trends of consolidation ; 
characteristics of layered systems under construction loading | 
obtained by mathematical analyses are within reasonable error 
as s>een in all the above experimental verifications. 
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Ins oantaneous loading 

Two and three layer systenis : A number of interesting trends 
and some fallacies in computing the time rates of settlement 
were demonstrated in Chapter III relating the effects of 
location of drainage boundary to the relative values of permeabi 
lity and compressibility coefficients in two and three- layer 
systems. Due to experimental limitations it has not been 
possible to demonstrate these trends. It would have been 
interesting to obtain settlements of individual layers even 
thou^ it will be nearly impossible to control the individual 
consolidation parameters of the layers. However, a comparison 
of the total system in experiment with the analytical results 
is useful to demonstrate the trends. 

The results of two and three layer systems under 
instantaneous loading are presented in Pig. 5.24 to 5.30. 

In all the cases -it is seen that the analytical solution 
obtained by numerical method compares more favourably than 
the transformation method with the experimental results. In 
some cases the transformation method yields very crude results 
(Pigs. 5.26, 5.27 and 5.28). 

CONCnJSIONS 

A simple experimental set-up is deviced for applying 
gradually increasing construction type of loading in one- 
dimenstonal consolidation test. Special oedometer rings are 
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used to conduct tests on layered systems. Using four different 
types of soils, one-dimensional consolidation tests on layered 
systems up to three layers arc conducted under instantaneous as 
well as construction loading. 

The general trends of consolidation characteristics 
of layered systems under instantaneous as well as construction 
loading predicted by analytical procedures as discussed in the 
earlier chapters are found to be within reasonable limites of 
errors when compared with the experimental results-. 

A simplified method of analysis of the layered system 
under instantaneous loading, using the method of transf ermation 
of individual layers, results in substantial errors. 
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CHAETER VI 


coNCiusioNs nmoMmmAiiom 

COl'JCIUSIONS 

An attempt is made in this research to obtain solutions 
to several complex field conditions and loading environments so 
as to arrive at certain general conclusions on the consolidation 
behaviour of soils subject to these conditions. A limited number 
of laboratory experiments duplicating some of the above situations 
have been conducted and the results are compared with the 
analytical solutions. Ihe following conclusions summarise 
the findings of this investigation: 

For a thick clay layer even though the coefficient of 
consolidation C. is coriStant throughout the depth, the rate of 
consolidation is veiy much affected by the individual variation 
of k and m with respect to depth. 

The location of the drainage boundary has significant 
Influence on iiie consolidation characteristics of thick clay 
layer with single face drainage. The consolidation is faster 
when the drainage boundary is located on the side of higher k 
and m values . 

A set of design charts for a thick clay layer under 
construction type of loading are presented. If appropriate 
field situations could be idealised under the various cases 
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pcrcscnted, the design charts will be useful in predicting the 
field settlement. 

A new formulation taking into account of the general 
variation of the consolidation parameters is presented. This 
forms a non-linear theory of a consolidating non-homogene ous 
thick clay layer with the parameters k and m varying during 
the process of consolidation. IThe solution is obtained by 
step-by-step numerical integration. 

The effect of decrease of k during consolidation is 
found to be more than that of m. 

Even though the coefficient of consolidation is constant 
throu^out the depth and during consolidatiorw the rate of 
consolidation is influenced by the individual variation of k 
and m. 


. In layered systems the effects of extreme ranges of 
parameters of individual layers are studied. It is noted that 
a total dependence on the value of C for evaluating the consoli- 
dation characteristics of layered system will result in erroneous 
estimations. The relative magnitudes of k and m have much 
effects on the eonsolidation behaviour of layered systems. 
Therefore a layered system v/ith layers of same coefficient of 
consolidation cannot be treated as single layer always. 
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A set of design charts in a two- layer system for a 
limited range of soil parameters is presented. This will be 
of use to predict the settlement rates in the field. 

Closed form solutions are presented for single and 
two laj'-er system under construction type of loading for single 
face as well as two face drainage. A three-layer system under 
construction loading is analysed by finite difference method. 

The solutions are illustrated with examples. 

A simple experimental set-up is deviced for testing 
layered systems under construction type of loading. Using 
four different ytypes of soils, one-dimensional consolidation 
tests on layered systems up to tlnree layers are conducted. 

The general trends of consolidation characteristics of layered 
systems under instantaneous as well as construction loading 
predicted by analytical procedures are found t o be within 
reasonable limits of errors when compared with the experimental 
results . ' 

In general, the simplified analysis yield very 
approximate estimates as compared with the mere refined analysis 
presented here. 


EEC OlffiEimAT I Oils 

The follov/ing topics are recommended fer further 
investigati on . 
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Study of the variation of the consolidation parameters 
assuming more generalised functions such as; 


k 





M 


m 




T/here N, M, R and S are constants . y 

Experimental study in field to correlate the actual 
variation of consolidation parameters. 

Experimental study on the consolidation behaviour of 
individual layers in layered systems v/ith pore pressure measure- 
ments reflecting the effect of the location of drainage bound- 
aries, the relative magnitudes of the consolidation parameters^ 
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APPEIIDIX 


CON^^RCtENCE OE the HUBtERICAL RESULTS 
OE THE BTON-LIHE.'iR EQUATION 2.30 

The numerical compu.tation of the non-linear equation 
2.30 for- the particular case discussed, the convergence of the 
results are obtained as follows: 

Eor : o4 = -0.80 

|3 = 4.00 

= 0.50 

0.25 

At time factor T = 0.01 


VALUE CE PORE PRESSURE RATIO 


At Depth 

X 

Eor time interval = .0001 

Eor time interval = .001 

0 

.000000 

.000000 

0.1 

.636742 

.645151 

0.2 

.925340 

.93517 

0.3 

.991553 

.994327 

0.4 

.999446 

.999755 

0.5 

.999977 

.999995 

0.6 

.999999 

1 .000000 

0.7 

.999988 

.999993 

0.8 

.999600 

.999703 

0.9 

.987463 

.988543 

1 .0 

.000000. 

.000000 
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Steps in solving two layer system 
under construction loading 


The equations to be solved are; 


2 

9 u 


1 



+ d. 



( 1 ) 



+ d. 



( 2 ) 


The boundary conditions are as given by Eq. 4.20. 

Using the Laplace Transformation defined as; 

oo 

= fe~^^ u^ dt = ^5) 

0 

OO 

= / U2 dt = V2 (4) 

0 

The Eqs . (1) and (2) aro transformed as; 


Lju2^ 


d v^ d. 

7^ + — = sv^ 
dp 


d Vp dp 


dp 


__ -j- US Yo 

2 s r * 2 


(5) 

( 6 ) 


How the transformed boundary conditions to be satisfied by the 
Eqs. (5) and (6) are; 



At p = 0 


(7) 
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P = 1 

p = 1 


^2 

= 0 

(8) 

^1 

= ^2 

(9) 

flv^ 

3P 

9v, 

= ko — ^ 

2 9P 

(10) 


Tile solutions of the ordinary differential equations 
(5) and (6) satisfying the boundary conditions (7) and (8) only 
are ; 


V, 


d^ < j 

= — 1 1 - Cos h p + M Sin h v/s p j (11 ) 


V2 = -^ U - 


Cos h (l+v--p) + N Sin h (1+V-p) y/s ^ 

( 12 ) 


V/here M and F are constants to be evaluated satisfying the 
b oundary c ondit i ons ( 9 ) and (10). 


Using the notation A = >/s and substituting the boundary conditions 
(9) and (10): 


-trSin h 1 Sin h Afiv+ Cos h A Cos h A^V-1 
Sinh A Cos h AyAV+ Cos h A Sin h Apv 

■ (13) 

< 3 - Cos h A Cos h A}iv- Sin h A Sin h ApV-d" 
Sin h A Cos h Apv+ 'Cos h A Sin h Ajjv 

( 14 ) 
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Therefore , 


= 




1 n-_i, -w , ~'*‘Sinh A Sinh Ahv+ Cosh Cosh At^v-I 
1-Cosh ^P+— SiHTA'-'Cosh Afiv+ C C^^;hT"Sinh'ljw 


Apj 


(15) 


To find out the final solution, the Inversion Theorem of the 
Laplace Transf orEBtion is used; 

I + i'" ' 

I _100 

■'iTaere is to he so large to include all the singularities of 
the integrand in Eq. (16). 

The line integral of Eq. (1b) is found out by using the well- 
knovm Residue Theorem. 


The integrand as per Eq. (14). has simple poles at s = 0, 

2 

and s = - A^ where A^ are the roots of the following equation 


Sin A Cos A^V +<5^003 A Sin A^V = o (17) 

Also the above equation has all roots real and simple (8, 9, 
10). Therefore: 


Res 


s=0 




i 


Vi (I) = 2^ X 2« 


(18) 
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To find the residues; 

The dinominator is obtained as ; 

ds(Sin A Cos A|i^)+(rCos A Sin AjjV) ( 19 ) 


= A jptSin A Sin A|..iV - Cos A Cos ApV 

-(5'pvCos A Cos Apv + Sin A Sin Aj 4 V|* 

= Apv ^ Sin A Sin A|A\) ~ (T Cos A Cos Ap'vJ 

- A ^ - <rsin A Sin Apv + Cos A Cos Ap V j 
Using Eq . ( 22 ) the above expression is simplified as follov/s 

. = Apv|sin A Sin A^V + | 

_ js I Sin^^O£s_AJll^ ^ 


_ AKV 
Sin A 


Sin A 


. Cos AHV 
CosA 


AftV 
Sin Apv 


Sin A + A 


Sin A^v 
Sin A 


= ^ Sin^A^v j (20) 

The Numerator is obtained as; 

The contribution by the first two terms of the 
Eq. (15) is zero, so the remaining terms are; 

^-crSin A Sin Apv + Cos A Cos Apv - l| Sin Ap 

( 21 ) 
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Using "le Eq.(16) and the expression (20), the Residue v/hen 
s = 0 i . obtained as ; 


0"Sin^ Apv + Sin A Sin AP'' a„ 

^ ^ — bin Ap 

A-" |.tVSin A +<?'Sin'^ Ajj.v 


( 22 ) 


2 

And the Residue at s = -A^. is obtained as; 


E —1 ' g"Sin‘^AHVi- Sin A Sin Apv c..- 

pvSia^A +<TSin^Apv) 


-A^T 

Sin Ap e (23) 


Therefore, the final solution (p, T) is put in the forn as 
given in Eq, 4.2'i . 
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